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ABSTRACT 

We investigate the chemo-dynamical effects of multiple supernova (SN) explosions in the central region of 
primordial galaxies (e.g. M ~ 1O 8 M halo at z ~ 10) using three-dimensional hydrodynamical simulations of the 
inhomogenous interstellar medium down to parsec-scales. We find that the final protogalactic structure and metal 
distribution depend strongly on the number of SNe. Specifically, 1) 1000 SNe after an instantaneous burst of star 
formation are sufficient to almost completely blow away the gas in these systems, whereas 2) 100 SN explosions 
trigger the collapse of the protogalactic cloud, leading to the formation of a cold, dense clumpy disk (n > 300 
cm -3 ) with metallicity, Z = 4 x 10~ 4 Z Q . These results imply that the metallicity of the "second generation" of stars 
could be Z ~ 1O~ 4 Z , and that the environment to form metal-free stars in protogalaxies may be lost relatively 
quickly (< 10 7 yr) after the first burst of Z = star formation. The recently discovered ultra-metal-poor star 
(Christlieb et al. 2002) may represent surviving members of such second-generation star formation. 

Subject headings: ISM: structure, kinematics and dynamics — ISM: abundances — cosmology: theory — 
galaxies: formation — methods: numerical 



1. INTRODUCTION 

The forefront problems of cosmology and star formation 
have become increasingly intertwined in the last decade. This 
has been driven especially by studies of the feedback from the 
first stars and supernovae in high-redshift (z) protogalaxies on 
the intergalactic medium (IGM). The collapse of primordial gas 
in dark matter (DM) halos followed by the formation of the first 
metal-free stars has been extensively studied by several groups 
(e.g. Omukai & Nishi 1998; Abel, Bryan, & Norman 2000; 
Bromm, Coppi, & Larson 2002). Examining the feedback from 
the first generation of supernovae (SNe) is essential for under- 
standing the metal enrichment and evolution of the interstellar 
medium (ISM), as well as the structure and subsequent star for- 
mation in the host protogalaxies. Some of the key physical is- 
sues in this problem are the propagation of metal-rich SN gas 
through the inhomogeneous ISM, and the subsequent mixing 
and diffusion of metals in the multi-phase ISM. Mori, Ferrara, 
& Madau (2002; hereafter MFM02), using three-dimensional 
hydrodynamical simulations, found that a significant fraction of 
the metal-enriched gas in a 1O 8 M halo at z = 9 is blown away, 
and that multiple SN explosions in primordial galaxies can play 
an important role in the metal pollution of the high-z IGM. 
However, the early stages of the interaction between multiple 
blast waves from SNe in an inhomogeneous ISM are important 
for the subsequent evolution of the ISM in protogalaxies. From 
this point of view, the minimum resolved scale in MFM02- 22 
pc- is not fine enough to resolve the interactions between the 
blast waves and the clumpy ISM. Furthermore, self-gravity of 
the gas and radiative cooling below 10 4 K, both of which are 
not considered in MFM02, are crucial factors in determining 
the internal structure of primordial galaxies after the initial SN 
explosions. This is especially relevant for cases in which the 
initial starburst is not strong enough to blow away the proto- 
galactic gas. 

In this paper, we report on the effects of an initial metal-free 
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starburst in primordial galaxies, using three-dimensional hy- 
drodynamical simulations with spatial resolution 3.9 pc, which 
take into account the self-gravity of the gas and radiative cool- 
ing between 20 K and 10 8 K. We focus here on the simplest 
scenario, in which all the SNe are coeval and spatially concen- 
trated; variations on these initial conditions will be considered 
in our future work. We explore feedback from relatively small 
numbers of SNe, in order to see what kind of ISM structure and 
metal distribution are achieved in primordial galaxies just after 
the^zraf SN explosions of first-generation stars. Studying SN 
feedback at these ISM scales has important consequences not 
only for the metal pollution of the high-z IGM, but also for the 
duration of metal-free star formation in the early universe. 

2. NUMERICAL METHODS AND MODELS 

The hydrodynamical schemes which we use are described in 
Wada & Norman (2001) and Wada (2001). We follow explic- 
itly the diffusion of metals ejected from SNe into the ISM by 
solving the equation similar to that of mass conservation and 
assuming that the metal component has the same local veloc- 
ity as the fluid. We calculate the gas cooling rate using cooling 
functions (Spaans, private communication) for 20 < T g < 10 8 K 
and for gas metallicities in the range Z = - 1Z Q . In the pri- 
mordial gas, H2 is the main coolant below 10 3 K. The H2 and 
HD abundances are determined through a chemical network. 
The diffusion of metals ejected from each SN into the ISM is 
explicitly followed. Therefore the gas metallicity is a function 
of position and time. The radiative cooling for the gas in each 
grid cell is determined for the local metallicity (Goldsmith & 
Langer 1978; Sutherland & Dopita 1993; Spaans 1996). Note 
that in reality the CMB sets a physical temperature floor for the 
radiative cooling, e.g. ~ 30 K at z ~ 10 (Bromm et al. 2002), 
although the gas temperature below ~ 100 K is not very impor- 
tant for the evolution of the blast waves. 

The hydrodynamic part of the basic equations is solved 
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by AUSM (Advection Upstream Splitting Method). We used 
512 2 x 128 equally spaced Cartesian grid points covering a vol- 
ume of 2 kpc x 2 kpc x 1 kpc. Therefore, the spatial resolution 
is 3.9 pc and 7.8 pc for the x—y and z directions respectively. 
The Poisson equation is solved to calculate the self-gravity of 
the gas using the FFT The second-order leap-frog method is 
used for the time integration. We adopt the implicit time inte- 
gration for the cooling term. 

To achieve density and temperature distribution in a dark 
matter potential for the initial condition, we performed a pre- 
calculation in which a geometrically thin, dense disk is evolved 
in the spherically symmetric, time-independent DM potential, 
$ DM . We assumed $ DM = -(27/4) 1/2 v 2 /(r 2 + a 2 ) 1 ' 2 , where 
a = 267 pc is the core radius of the potential, and v c = 1 1 km s" 1 
is the maximum rotational velocity. The mass of the dark mat- 
ter in the central part of the protogalaxy is therefore 2 x 10 7 M Q . 
The total mass and initial temperature of the gas are 1O 7 M and 
20 K. The gas is uniformly distributed in a thin axisymmetric 
disk, whose radius and scale height are respectively 400 pc and 
8 pc. The total DM mass of ~ 10 8 M Q corresponds to a virial 
radius of ~ 1 kpc at z ~ 10 (MFM02; Madau, Ferrara, & Rees 
2001). Note that the very first stars could form in smaller sys- 
tems of mass ~ 10 5_6 M Q at higher redshift (z^ 20-30) as stud- 
ied in Bromm, Coppi, & Larson (2002). Moreover, in the low- 
mass galaxies considered here, star formation followed by SNe 
might already have occurred during the initial collapse of the 
DM component. For simplicity, and owing to the short dynam- 
ical time scale of the gas component after the SN explosions, we 
assume here that the DM potential is time-independent. A more 
self-consistent treatment will be considered in future work. 

In the upper panels of Fig. 1 (t = 0), the density and tempera- 
ture distribution of the gas after virialization are shown. During 
the virialization process, the density distribution becomes al- 
most spherical, although some dense clumps are formed in the 
central region where the temperature is less than 10 3 K. The 
mean density and temperature within 100 pc from the center 
are 3 cm" 3 and 100 K. The total gas mass in the computing box 
is 5.6 x 10 6 M Q . Gas in the cold phase (T g < 100 K) occupies 
about 1.4% of the volume, and has a mass of 1.9 x 10 6 M Q . The 
maximum density in the dense clumps is 20 cm" 3 (Jeans length 
~ 18 pc). We use this density, temperature, and velocity distri- 
bution as the initial condition. 

In order to assess the chemo-dynamical effects of the first 
SNe, some assumptions must be made about the initial mass 
function (IMF) of the first stars. There are significant uncer- 
tainties associated with this, since the processes that determine 
the stellar IMF, even in the local universe, are at present unclear. 
Some recent theoretical studies (e.g., Abel, Bryan, & Norman 
2000; Bromm, Coppi, & Larson 2002) indicate that the primor- 
dial IMF may have been weighted towards masses exceeding 
100 M Q . If this were true, the metal yield per SN could poten- 
tially be much higher than that for the "usual" Type II SN. How- 
ever, the duration and universality of the processes that lead to 
such an IMF are currently unknown. Given the unresolved na- 
ture of this topic, we assume a present-day Salpeter IMF over a 
stellar mass range of 1 - 100M Q in this work. We include metal 
generation only from Z = stars of mass 8 - 40M Q (Woosley 
& Weaver 1995), assuming that stars more massive than 40 
Mq implode directly into a black hole without mass/metal ejec- 
tion into the ISM. For these assumptions, each unit mass of gas 
that forms stars produces total IMF-integrated masses in metals 
and stellar ejecta of 0.00663 and 0.1863 after the SN explosion. 



Furthermore, for the above assumptions, there is one SN for ev- 
ery ~ 100 Mq of stars formed. Thus, the total metal yield and 
ejecta mass per SN on average are 0.663 M Q and 18.63 Mq. 
We input the energy (10 51 erg per SN) as thermal energy and 
mass/metal yields from the SNe at time t = 0. We also explored 
cases in which the SNe are not coeval, and occur over a period 
of <~ 1 Myr. However, no essential differences from the models 
discussed in this paper were found. 

Since the IMF and star formation efficiency in primordial 
galaxies are poorly understood, we take the number of SNe, 
7Vsn> as a f ree parameter, and we explore two models here: 
model A (7V SN = 10 3 ) and model B (7V SN = 10 2 ). We put in N m 
SNe randomly within the central 100 pc, since some of the ini- 
tial star formation is expected in the dense, cold clumps seen 
in Fig. 1 (f = 0). The energy and ejected masses in metals and 
gas from each SN are injected into a randomly selected grid 
cell, around which a blast wave is then expanded. This wave is 
not necessarily spherical, owing to the inhomogeneous density 
field and interactions with other SNe. 

A priori, we expect that if the total input energy is sufficiently 
larger than the binding energy, then the blast waves will blow 
the gas away, as has been predicted in many papers (e.g. Larson 
1974, and references above). On the other hand, if the total SN 
energy is not efficiently converted into the ISM kinetic energy, 
the multiple blast waves would form a cavity in the protogalac- 
tic cloud. We would then expect the whole system to become 
dynamically unstable and collapse. This can be demonstrated 
in an analytic estimate as follows. Using a similarity solution 
of a blast wave in a radial density profile going as, p oc r" 2 , 

the maximum radius of blast waves, r s = (6E / p) l / 5 T^ 5 for 
7 = 5/3 (Ostriker & McKee 1988), where E is the input en- 
ergy, p is the average density, and t c is the cooling time of 
the hot gas in the cavity generated by the SNe. In our work 
here, t c ~ 2 Myr, so that the cavity radius would be ~ 200 
pc (£/O.l) 1/5 (7V SN /lOO) 1/5 (p/O.O3M pc- 3 )" 1 / 5 (T c /2Myr) 2 / 5 , 
where E = e/VsnIO 51 erg, and e is the heating efficiency (the 
fraction of the total SN energy that is converted to the kinetic 
energy of the blast waves). When the hot gas in the cavity 
cools, we expect the protogalactic gas outside the cavity to lose 
pressure support and to begin falling towards the galactic cen- 
ter. The gas mass inside the cavity (r < r s = 200pc) is about 
10 6 M Q . If all the gas in the protogalaxy, about 5.6 x 1O 6 M , 
were to collapse and mix uniformly with the hot metal-enriched 
cavity, then the resulting average gas metallicity would be about 
100 SN x (0.7 M /SN)/(6.6 x 1O 6 M ) - 5 x 10"%. 

In reality, the above estimates will be affected by the evolu- 
tion of the multiple blast waves in an inhomogeneous medium, 
the fraction of protogalactic gas that eventually collapses, and 
by the heating efficiency, e, which in this situation is unknown. 
In addition, the mixing of metals in the collapsing protogalaxy 
will be nonuniform. In our numerical method here, we do not 
have to make any specific assumptions for these factors. In- 
stead, we explicitly solve for the metal diffusion, and the in- 
teraction between the multiple blast waves and the inhomoge- 
neous ISM. 

3. RESULTS 

The evolution of the density and temperature distributions for 
model A are shown in Figure 1. The multiple SN explosions 
evolve as one big blastwave, which reaches the outer boundary 
of the z-direction at t = 2.1 Myr, and those in the x— and y— 
directions at t ~ 6 Myr. As seen in the figures, the explosions 
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are roughly spherical, but dense 'spikes' are formed, caused 
by non-linear instabilities in the decelerated shocks (Yoshida 
& Habe 1992). The hot gases in the cavity quickly cool from 
- 10 7 K at t = 2.1 Myr to - 10 3 -10 6 K at t = 5.8 Myr. At t = 3 
Myr, the gas velocity is 2-4 times larger than the escape velocity 
in the whole region. This implies that 1000 SNe are sufficient 
to blow away the gas in the protogalaxy in <~ 10 7 Myr. 

On the other hand, in model B, the protogalaxy survives the 
initial starburst. The 100 SNe generate a cavity in the central re- 
gion of the protogalaxy, but it does not exceed sizes of r <~ 300 
pc. After t ~ 20 Myr, the cavity collapses, which in turn trig- 
gers the collapse of the protogalactic gas cloud. Eventually a 
clumpy disk-like structure, whose radius is about 200 pc, is 
formed in the center (Fig. 2). The y — z cross section reveals 
a 'chimney' -like structure. The diffuse gases seen in Fig. 2 are 
still falling in towards the central dense disk. At t = 52.6 Myr, 
the mass of dense gas (n > 30 and 300 cm" 3 ) is 1.6xl0 6 M© 
and 7.2xlO 5 M , which are 36% and 16% respectively of the 
total gas mass in the box at that time. The maximum density is 
3.4xl0 4 cm" 3 . 

Fig. 2 (c) shows that the metallicity distribution in the 'chim- 
ney' is about (10~ 3 -10~ 2 ) Z , and that the diffuse halo is almost 
metal-free. The average metallicities in the very dense regions 
of n > 30 and 300 cm" 3 are 6.5x10"% and 4.1xlO- 4 Z . 
These are in contrast to Z > 1O" 2 Z in the relatively low density 
regions between the dense clumps in the central disky structure, 
owing to the formation of these clumps through the accretion of 
mainly metal-poor or metal-free gas. 

As one may expect, the two models have different conse- 
quences for the global ejection of metals. Fig. 3 shows the 
evolution of the radial metal mass distribution with time in the 
two models. Almost all the metals (6 x 1O 2 M ) in model A are 
ejected out of the central protogalactic region, whereas most of 
the metals in model B return to r < 200 pc in 50 Myr. As a re- 
sult, the central region in model B is ten times more metal-rich 
than in model A. 

In addition, the efficiency of converting the SN energy to the 
kinetic energy of the ISM is time-dependent. In model B, this 
efficiency, defined as the ratio of the total ISM kinetic energy 
to total input energy, is 0.12 at t = 1 Myr, and 7 x 10" 3 at t = 10 
Myr. 

In summary, our results in the context of the simplified mod- 
els we have considered here- SN feedback from an instanta- 
neous starburst in protogalaxies of mass about 1O 7 -1O 8 M - 
are as follows. First, protogalaxies which undergo bursts with 
more than 10 3 SNe are not likely to retain the generated met- 
als. The ultimate fate of the ejected gas and metals in this case 
is unclear, as they have not been tracked for sufficiently long 
timescales here. They may eventually escape to the IGM, or 
fall back partially to the galactic center, forming a dense core 
as found in the large dynamic range simulations by MFM02. 
Second, we find that less active starbursts (TVsn ~ 100) could 
lead to faster metal re-incorporation in the gas, from which 
"second generation" stars with Z <~ 1O~ 4 Z could potentially 
form. The value of this transition metallicity is indicative rather 
than exact, given that we do not explicitly track the process of 
star formation itself, and that mixing processes depend on in- 
stabilities that may not be entirely resolved by our simulations. 
Nevertheless, this metallicity is interesting in terms of the re- 
cent discovery of an ultra-metal-poor star with [Fe/H]= -5.3 
(Christlieb et al. 2002). The abundance pattern of elements in 
this object suggests that it may have formed from a gas cloud 



that was pre-enriched by Type II SNe. Such a star could rep- 
resent the surviving members of the second-generation stellar 
population in our scenario's model B, but with ten times fewer 
initial SNe. In addition, our derived value for the transition 
metallicity is in agreement with the numerical estimates in re- 
cent works (Bromm et al. 2001; Mackey, Bromm, & Hernquist 
2003; Schneider, et al. 2002), for the critical level of enrich- 
ment that must occur in order for star formation in primordial 
gas to evolve from a hypothesized top-heavy stellar IMF to a 
present-day IMF. Future observations of ultra-low-Z stars will 
strongly test these theoretical predictions for the metallicity and 
masses of second generation stars. 

4. DISCUSSION 

In this paper, we have begun an investigation of SN feedback 
in high-z protogalaxies. Our methods contain significant im- 
provements over recent works, such as the inclusion of gas self- 
gravity in 3D and improved spatial resolution. These are nec- 
essary to quantitatively understand the effects of multiple SNe 
from the first stars at scales which are important for the ISM, 
which influences the cosmological relevance of early SNe for 
the IGM. We find that the fate of the generated metals- whether 
they are transported as a hot gas phase to the IGM or reincor- 
porated quickly into cold starforming gas in the protogalaxy- 
depends critically on the number of initial SNe. This has broad 
implications for the transition metallicity at early epochs from 
metal-free to "second generation" star formation, and the rapid- 
ity of this process. The duration of metal-free star formation is 
of critical importance to theoretical studies of the H I and He II 
reionization of the IGM (see, e.g., Venkatesan, Tumlinson, & 
Shull 2003), as well as for observational searches for such stars 
in the Galactic halo or at high redshift (Tumlinson, Shull, & 
Venkatesan 2003, and references therein). 

In the initial condition, the total mass of cold (T g < 100 K), 
dense gas (n > 3 cm" 3 ) is 3 x 1O 5 M . The gas mass converted 
to stars in model A, with 10 3 SNe, is about 1O 5 M . Therefore, 
if the star formation efficiency is larger than ~ 0.3, the first ob- 
jects formed in the protogalaxies are almost pure stellar systems 
of mass ~ 1O 5 M , and the remaining baryonic mass is ejected 
from the galaxy. These metal-free stellar systems would not be 
bound, because a large fraction of the initial gas mass in the 
protogalaxy's central region is lost due to strong winds from 
the starburst's SNe. On the other hand, if the initial star forma- 
tion efficiency is of order a few percent, the first starburst does 
not inhibit the subsequent creation of new star formation sites. 
If the star formation efficiency remains low, the metals from 
SNe cannot be expelled from the protogalaxies. In this case, 
metal-enriched stellar clusters of mass ~ 10 5 - 1O 6 M even- 
tually form, accompanied by a diffuse gaseous halo. These 
systems can be the building blocks of larger galaxies, or the 
progenitors of globular clusters. In other words, suppose that 
the star formation efficiency is constant, e.g. 0.1, among pro- 
togalaxies. The smaller systems (M tota i < 1O 7 M ) would then 
dominate the metal enrichment of the IGM, which is consistent 
with the findings of Mac Low & Ferrara (1999), and Ricotti, 
Gnedin & Shull (2002). The fate of the metals generated by the 
first stars is clearly linked closely to the mass and star formation 
histories of the host protogalaxies. In this work, we have con- 
sidered a simple burst scenario in which the SNe were coeval 
and spatially concentrated, in order to focus on the ISM gas 
dynamics. We will examine the importance of varying these 
and other factors, such as the stellar IMF, in more detail in the 
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future. 
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FIG. 1. — Log-scaled density (Mq pc 3 ) and temperature (K) distributions of the initial condition (f = 0), and at two snapshots (t =2.1 Myr and 5.8 Myr) of model 
A (Msn = 10 3 ). (a) The x—y cross sections of density, (b) the y-z cross sections of density, and (c) same as (b), but for temperature. 



FIG. 2. — Same as Figure 1, but for model B (Wsn = 10 2 ) at t = 56.2 Myr. The lower right panel shows the y-z cross section of the metallicity distribution in the 
central region. The plotted range of metallicities is Z = lCr 4 Z0-lfr 1 Z0. 
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